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We derive a equation of motion for the velocity fluctuations 
of a 2D deformed trapped Bose gas just above the critical tem- 
perature in the hydrodynamical regime. From this equation, 
we calculate the eigenfrequencies and the corresponding den- 
sity fluctuations for a few low- lying excitation modes. Using 
the method of averages, we derive a dispersion relation in a 
deformed trap at very high temperature that interpolates be- 
tween the coUisionless and hydrodynamic regimes. We make 
use of this dispersion relation to calculate the frequencies and 
the damping rates for monopole and quadrupole mode in both 
the regimes. We also discuss the time evolution of the wave 
packet width of a Bose gas in a time independent as well as 
time dependent trap. 

PACS numbers:03.75.Fi,05.30.Jp,67.40.Db 



I. INTRODUCTION 

There has been renewed interest in the Bose-Einstein 
condensation(BEC) after its experimental demonstration 
by Anderson et.al |l) in a magnetically trapped rubidium 
gas. In other experiments the Bose-Einstein condensate 
of trapped lithium and sodium vapor is also ob- 
served. There has been much interest in the theoretical 
understanding of this system Q . The low- lying collective 
excitations above the critical temperature in the hydro- 
dynamic region has been discussed by Griffin et.al H] 
using the kinetic theory. Damping of the hydrodynamic 
modes in a trapped Bose gas above the Bose-Einstein 
critical temperature (Tc) is also discussed by Odelin et.al 
and Kavoulakis et.al |Q. Frequencies of the low- lying 
excitation modes at T = have been discussed by S. 
Stringari ||]. 

After the discovery of BEG in a trapped alkali atom, 
the influences of the dimension of a Bose systems has 
been a subject of extensive studies 1^. In our present 
technology one can frezze the motion of the trapped par- 
ticles in one direction to create a quasi-2D Bose gas. In 
the frozen direction the particles executes the zero point 
motion. To achive this quasi-2D system, the frequency in 
the frozen direction should be much larger than the fre- 
quency in the X-Y plane and the mean field interactions 
among the particle. It has been shown by V. Baganato 
et.al that for an ideal 2D Bose gas under harmonic 
trap a macroscopic ocupation of the ground state can 
exist at temperature T < Tc = xf^^- A = ^ is 



the deformation parameter. Some experiment shows the 
possibility of creating a quasi-2D trapped Bose gas [ pH . 

It is well known that the excitation frequencies for 
monopole and quadrupole modes are 2ujo and v^^o re- 
spectively in a 2D isotropic trapped Bose gas. Using the 
approximation, Uz » uj, the dispersion relation of the 
excitation frequencies does not produce the cor- 

rect frequencies for monopole and quadrupole modes in 
a 2D trapped Bose system. There has been no system- 
atic study on the collective excitations of a 2D deformed 
trapped Bose gas above the critical temperature. 

The purpose of this paper is to give analytic results 
for the dispersion law of low-lying collective modes in 2D 
deformed trapped Bose gas and their damping rates in 
both regime, hydrodynamic and coUisionless. 

Above the critical temperature, one can distinguish 
two regimes, the hydrodynamic (coUisional) one where 
collisions ensure the local thermal equlibrium and col- 
lisionless where the motion is described by the single 
particle hamiltonian. In the hydrodynamic region, the 
characteristic mode frequency is small compared to the 
collision frequecy {lut << 1). In the coUisionless region 
{ujT >> 1 ), the collision are not important. 

The paper is organised as follows. We derive in sec[II] 
a closed equation of motion for the velocity fluctuations 
of a 2D deformed trapped Bose gas just above the crit- 
ical temperature (T > Tc) using the kinetic theory. We 
make use of this equation in sec [III] to calculate the ex- 
citation frequencies for a few low-lying collective modes 
and the corresponding density fluctuations. In sec [IV], 
we derive a dispersion relation in a 2D deformed trap 
at very high temperature using the method of averages 
that interpolates between the coUisionless and hydrody- 
namic regimes. From this dispersion relation, we calcu- 
late the eigenfrequencies and damping rates for monopole 
and quadrupole mode. We discuss the evolution of the 
wave packet width of a Bose gas in a time independent 
as well as time dependent trap. In sec[V], we presents a 
summary of our work. 



II. HYDRODYNAMIC EQUATION OF MOTION 
FOR THE VELOCITY FLUCTUATIONS 

We shall discuss the collective modes of a 2D deformed 
trapped Bose gas in the hydrodynamic regime just above 
the critical temperature T > Tc using the kinetic theory. 
In the low energy excitations, we can use the semiclassical 
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approximation for the dynamics of a Bose gas, using the 
following Boltzmann equation for the phase space 
distribution function f{r,p,t) 



dl 

at 



(1) 



where IcoU is the collisional integral and F ~ —WUo(r}. 
The trap potential is Uo{r) = im(ti;^a;^ + w^y^). 

In the hydrodynamic regime, collisions ensures the lo- 
cal thermodynamic equilibrium. To lowest order, the per- 
turbed distribution function produced by a slowly vary- 
ing external potential is the equlibrium Bose distribution 
function 



f{r,p,t)^[exp{f3{r,t)rj{r,t)-l]- 



V{r, t) = - /i(r, t) 



(2) 



(3) 



/i(r, t) is the chemical potential. The inverse temperature 
is [3{r,t) = Y^^^pj^- We are only interested in a small 
perturbations around the equilibrium states. 
The conservation laws are : 



i9n(r, t) 
dt 



dSv 



+ V.[noir)Svir,t)] = 



(4) 



mno(f)— = ~[VP(f,i) + n(f,t)Vi7o(r)] (5) 



dE{f, t) 
dt 



-V.[(P(f) + E{r))6v] ~ na5v.VUo{r) (6) 



These conservation laws are obtained from Eq. p) multi- 

-12 

plying by 1 , p, ^ and integrating the resulting equation 
over p. During collisions, the total number of particles 

-2 

N, momentum p, and energy ^ are conserved, so the 
collisional term vanishes. 5v is the velocity fluctuation 
around the equilibrium states. 

Using the quantum statistical mechanics, pressure and 
density can be written as 



P 



92{z) 
A2 



A2 



(7) 



(8) 



where gn{z) — ^^'^ Bose-Einstein functions. 

z{r,t) = e'^f'"'*)^''''*) is the local thermodynamic fugacity 
which is always less than one. A = 
de-Broglie wave length. 

One can easily get the relation. 



mkbT 



is the thermal 



P{r,t) = E{r,t) (9) 
in 2D . Using Eq. (^, Eq. can be written as 

dP{r,t) ^ _2v.[(Po(f)5z;(f,t)] - noSvif,t).VUoir) 



dt 



(10) 



Taking the time derivative of Eq. (|^) and using Eqs. 
(|) and (|l^), we get 

d'^Sv Pnir) 

= 2^^V[V.5v] ~ [V.<5w]VC/o(r) - V[<5w.V?7o(r)] 
at'^ n^yr) 

(11) 

The closed equation of motion for the velocity fluctua- 
tions has been derived by Griffin et.al jsj for 3D trapped 



Bose system. The term of (11) is associated with 

the Bose statistics. 

Without any external potential Uq — 0, the Eq. ([ll| ) 
becomes 



^ 2— 4d-V[V.(5w] 



dt^ 



m{r) 



(12) 



It has the plane wave solution with the dispersion relation 



2 - c^fc2. The sound velocity is - „2 



r=r or c 



A'o('-) 



where zq = e'^ and /io(r) ^fji-Uo{r). At 
high temperature (z < < 1 ) , the sound velocity becomes 
^2 _ afcsTh^ This sound velocity exactly matches with 
known result. 

From the continuity Eq. (0) we have. 



d5n{r, t) 
dt 



-(V.(5u)no(f) - 6v{f,t).Vno{r,t) (13) 



The density fluctuation is given by Sn{'r, t) — Sn{f)e 

In classical limit, the static density profile is rio(r) = 

„ / 2„2| 2 2\ 

no(r = 0)e 55 ^ where Q — kgT. 



III. EIGENFREQUENCIES AND THE 
CORRESPONDING DENSITY FLUCTUATIONS 
IN THE HYDRODYNAMIC REGIME 



1) The normal mode solution of (11) is 6v{r) — V(z'), 
here z = {x + iy) and I > 0. The excitation frequencies 
and the associated density fluctuations are = luj"^ , 
6n.j: Ci;2xz('"^)no(r) and u'^ ^ loI + {I - l)ujy , Suy 

For isotropic trap, the frequency is w = uJoVl- The 
corresponding density fluctuation is 5n{'r) ~ no(r)r'e''^. 
At r = there is no density fluctuation. There is a 
maximum density fluctuation at r = , / — ^ . 

2) The other solution of Eq. ( pi] ) is 6v{r) — V[ax'^ ± 
/3y2]. The positive sign is for the monopole mode and 
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the negative sign is for quadrupole mode. In a deformed 
trap, the excitation frequencies are 

- + ^l) ± V'9(c.2+,.2)2_32^2^] (14) 

For an isotropic tarp, it becomes lo = 2a;o or w = \/2wo- 
Hence in the anisotropic trap the monopole mode is cou- 
pled to the quadrupole mode. \i << Wj,, the low- 



est excitation frequency is = 




the lowest excitation frequency is uj — y ^ujy . The 

density fluctuation for the monopole mode is Sn[r) ^ 

[2 — ^ ]no{r) where as the density fluctuation 

for quadrupole mode is Sn{f) ~ (w^j/^ — Lu'^x'^)no{r). 

3) There is another quadrupole mode which has ve- 
locity field 6v{r) — V{xy). This is also called scissors 
mode [|l6|. The excitation frequency is uj"^ — uj^ ^ Uy 
and the corresponding density fluctuation is 6n(r) ^ 
(uj'^ + ujy)xyno{f). In isotropic trap lu'^ = 2uJq, which 
agrees with that for the scissors mode in hydrodynamic 
regime above |l|]. 



X2 = y^~ (19) 

X3 = xv^ + yvy (20) 

X4 = yVy - XVx (21) 

X5^vl+ vl (22) 

X6 = - vl (23) 

Using the Boltzmann kinetic equation (|l5|), we get the 
following closed set of equations. 

< Xi 2 < X5 > < Xi > +e < X2 > (24) 

< X2 >= 2 < X6 > < X2 > +e < Xi > (25) 

< X3 >= 2e < X4 > -2i < X3 > (26) 

<X4>=2e<X3>-2i<X4>-^^ (27) 



IV. METHOD OF AVERAGES 

At very high temperature, the dynamical behaviour 
of a dilute gas is described by the Boltzmann tarnsport 
equation. Here we include the collisional term in the 
Boltzmann transport equation and study the eigenfre- 
quencies for monopole and quadrupole mode using the 
method of averages Q . These two modes are coupled in 
a deformed trap. 

From Eq. (|l|), one can get the equations for the average 
of a dynamical quantity x(^j^) is [|l3|, 

< U.V.PX > - < — -VffX >=< ^coUX. > (15) 

at m 

where the average is taken in phase space and < x > can 
be written as 

< X >= J d^r(fvf[f, V, t)x{r, v) (16) 

< xlcoii > can be defined as 

< xlcoii >^ ^ J d^rtfvlxi + X2 - xi - X2]Icoii{f) 

(17) 

If X = + b{r).v + c{r)'ip , for elastic collision the 



collisional term is zero [|6|, |13|. a, b, and c are arbitrary 
functions of the position. 

Now we define the following quantities, 



2 I 2 

Xi = x +y 



(18) 



< X5 > = e < X6 > < X5 > -ri < X2 > (28) 
e2- 



< X6 > 



< Xi > 



£2+^2 

-et < Xi > -\ 7. — < X2 > 



e < X5 > - 



2 

< X6 > 



(29) 



- i < X6 > 



where double dot indicates the double derivative with 
respect to time, t = loI + lo^ and e — uil — uSy. The 
X6 is not a conserved quantity, so the collisional contri- 
bution comes only through the xe term. We have used 
the fact that < XqIcoU >= where r is the relax- 

ation time. This relaxation time r can be computed by 
a gaussian ansatz for the distribution function. The re- 
laxation time T is order of the inver se of the collision 
rate ^coii ^ n{0)vthcro, where Vth = \f^^- is the mean 

thermal velocity and n(0) — ^^^•^^ is the central den- 
sity for a quasi-2D system, is the osscilator length in 

4a^ 



the z-direction. Hence r 



ao = 87ra2 

is the 3-D scattering cross-section. It can be written in 
terms of Tc as 



1 



(30) 



The relaxation time r varies as vT in a quasi-2D where 
as in 3D it varies as TM . Now we are looking for a so- 
" ~ ' ■ " ' t ^ \Ye have the following 



lutions of Eqs. (24)-(p9|) as 
dispersion relation 
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LOT 



(31) 



Suliol] = 



This dispersion relation interpolates between the coUi- 
sionless and hydrodynamic regimes. In the hydrodynam- 
ical regime ( lot ^ ), the first term does not contribute. 

It gives W2 = 1 [3(^2 ^ ^2) ^ ^^9(6^2 + ^2)2 _ 32^2^2], 

This eigen frequencies exactly matches with Eq. (|l^), a 
result we found using the equation of motion for the ve- 
locity fluctuations even in deformed trap also. We have 
considered a few low energy excitation modes for which 
W.Sv is constant. The first term of the right-hand side 
of Eq. (11) does not contribute in the excitation spec- 
trum. Thats why the frequencies of these normal modes 
are same for a Bose gas just above Tc and at very high 
temperature. In pure coUisionless regime ( lot — > oo ), it 
gives LOc = 2lOx and ujc = 2a; j,. 

We can write phenomenological interpolation formula 
for the frequency and the damping rate of the modes in 
the following form d)- 0, 



UJ 



1 



(32) 



The imaginary part of the above equation gives for the 
damping rate 



(33) 



2 1 + {uotY 



where d = (w^ — uifj). In the hydrodynamic limit {ujt 
0), the damping rate is 



where as in the coUisionless region {lot ^ oo), 



r, 



CL 



2L0i 



(34) 



(35) 



The damping rate depends on the difference between the 
square of the frequencies in the coUisional and hydrody- 
namical regime. The damping rates can be calculated 
for different values of temperature, number of trapped 
atoms as well as of the trapping parameters and scat- 
tering length through the relaxation time r (^0|) . For 
monopole mode in an isotropic trap the difference d is 
zero. So there is no damping in the monopole mode in a 
2D isotropic trapped Bose system when the tcmpareture 
is very high. It was first shown by Boltzmann [l^ ] and 
later Odelin et.al Q in a 3D trapped Bose system at very 
high temperature. 

For isotropic harmonic trap, Eqs. (^^ - ( p9| ) decouples 
into two subsystem, one is for monopole and other one 
for quadrupole mode. The closed set of equations for 
monopole mode are 



< xi >= 2 < X5 > -2a;2 < xi > 



< X5 >= '^^Q <Xi> +2w^ < X3 > 



(36) 



(37) 



(38) 



There is no coUisional term in the above equations. So 
there is no damping for the monopole mode of a classical 
dilute gas confined in a isotropic trap. We are looking 



for solutions of Eqs 
The Eqs. (|3~ 



6D 




|) as e 



we get uo ~ 2uji 



'0- 



< Xi > 



1 



2<Xi > 



) can be re-written as 
(<Xi >') + 2uolxi - 



Xi 



(39) 



where Q = 2(< xi >< X5 > ^ < Xa is invariant 
quantity under time evolution. Define Xit) = ^< xi > 
which is the wave packet width and substituting it into 
Eq. !§) gives. 



X 



(40) 



This is a nonlinear singular Hill equation. The same 
equation is obtained at T = in 2D by Garcia RipoU 
et.al [|l5|. At equhbrium, = We linearized the 

Eq. (^ ) around the equilibrium point Xq, we get 



5X 



Auj16X = 



(41) 



One obtains an oscillation frequancy of the gas is w = 
2lo{), corresponding to the frequency of a single particle 
excitation in the gas. 

For time dependent trap, the equation of motion for 
the width of the wave packet is 



X + uol{t)X = ^ 



Q_ 

X^ 



(42) 



The general solution Q is X{t) = ^u^{t) + -§^v'^{t) 
where u{t) and v{t) are two linearly independent so- 
lutions of the equation p + LOQ{t)p = which satisfy 
u{to) = X{tQ), u{to) = X'(io), vito) = 0, v'ito) + 0. 
W is the Wronskian.This time dependent Hill equation 
( ^ ) can be solved explicitly only for a paricular choices 
of wo(i)- 

The closed set of equation for quadrupole mode in 
isotropic trap are 



< X2 > = 2 < X6 > -2c^o < X2 > 



< X4 > = -4cJo < X4 > - 



< X6 > 



(43) 



(44) 



< X6 > = 2^0^ < X2 > 



2uol < X6 > (45) 
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Solving this set of equation, we get damped quadrupole 
mode, 



LOT 







(46) 



In the hydrodynamic regime, the osciUation frequency is 
uj'^ = 2u!q where as in the cohisionless region, the fre- 
quency is just a single particle oscillator frequency. The 
damping rate can be calculated by using the Eqs. ( |3^ ) 
and @. 



V. SUMMARY 

In this work, we derived the equations of motion for 
velocity fluctuations of a Bose gas in a 2D deformed 
trap potential just above the critical temperature. When 
Uq = 0, it becomes a wave equation, from which we found 
the exact sound velocity at high temperature. We have 
also computed the frequency of the scissors mode in hy- 
drodynamic regime above Tc which agrees with the result 
obtained by Odelin et.al We have also calculated 

the frequencies for monopole and quadrupole mode and 
the corresponding density fluctuations in a deformed trap 
above T^. 

Using the method of averages, we obtained a disper- 
sion relation that interpolates between the cohisionless 
and hydrodynamic regimes at very high temperature. In 
a deformed trap as well as an isotropic trap, we have 
found the frequencies and the damping rates (in terms of 
the relaxation time) for monopole and quadrupole mode 
in both the regimes. In hydrodynamical regime, the ex- 
citation frequencies for monopole and quadrupole mode 
exactly matches with the previous result that we have 
found from equatin (|ll|). We have also shown that the 
relaxation time r varies as Vt in quasi-2D Bose gas. 

We have shown that there is no damping for monopole 
mode in a 2D isotropic trapped Bose gas when the tem- 
perature is very high. It was shown by Boltzmann [ p^ , 
later Odelin et.al |^ for 3D isotropic trap. 

We have discussed about the time evolution of the wave 
packet width of a Bose gas in a time independent as well 
as time dependent isotropic trap. It can be described by 
the solution of the Hill equation. 
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